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ABSTRACT 


The dynamic equations for a Liquid Metal Fast Breeder 
Nuclear Reactor (LMFBR) Power Plant are derived. These non- 
linear differential equations with typical design values for 
PmieesyStem Parameters are solved on a COMCOR Ci 5000 analog 
Computer for wie characteristic transients of the basic plant 
and the results compared with a recently proposed LMFBR design. 
A digital computer automatic control program is developed 
which uses, in addition to the ideas of derivative and pro- 
portional control, precomputed control trajectories. The 


digital controller design is implemented by applying FORTRAN 


TY? me tae 4. Vv 7 nes = 
Puno. Cauamnlnig co tne xbo 9SCC cigital computer Maemoe Toc 


tiveness of the strategy for the automatic controller is 
demonstrated by using various step changes in the load 
Semand . 

A special software routine was written to automatically 
scale the analog computer and proved to be quite uSeful in’ 
improving the operating efficiency of the hybrid computer 


samulation. 
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I. INTRODUCTION 


In view of the enormous increase Of world power consumption 
and the fast decrease of the resources which may be used to 
satisfy this growing power demand, it became necessary to 
develop and build nuclear power plants. The reactors used to 
generated the power for these plants are up to now almost 
exclusively thermal reactors using uranium as fisSionable 
material. There are 3 fissionable isotopes: U (235), Pu (239), 
me Ue ss)2 U (235) is found in natural uranium and can be 
used in this state or in an enriched form, i.e., its percentage 
1S increased beyond that in natural uranium. Pu (239) is 
meer eeturea Dy mammelear react1en from U (238). U (233) is 
Similarily manufactored from Th (232), an isotope which is also 
found in natural uranium. 

The worlds supply of nuclear fuel consists of 25x10° tous 
Geiaeurat uranium, mostly U (238), and 1x10° EGMemor COT Lunn, 
mostiy Th (232) [10]. The natural uranium consists of 99.298 
Weer) “mmemo./1% U (235). Since U (235) is the only naturally 
fissionable isotope the use of nuclear fuels would be limited 
wor0./1% of 25x10° tons. This means that, using only the 
present thermal type reactors with their inefficient nuclear 
economy a shortage of nuclear fuel is lightly to occur by the 
year 2000. 

Am increased muclear economy in chet un redctOr can we 


achieved by using fuel breeding. In this process the nonfission 


7 ae ee a 





capture of neutrons by a fertile U (238) nucleus results in 

the generation of fissionable Pu (239). Fertile Th (232) can 
also be converted to fissionable U (233). If more fissionable 
nuclei are produced than are consumed the reactor in which this 
process takes place is called a Breeder Reactor. If the fission 
process depends on fast neutrons as in the above described 
processes the reactor is a Fast Breeder Reactor (FBR). 

The control Problem of a FBR is particularly acute because 
of two important facts. First of all due to the high energies 
required of the neutrons in order to produce fission, the power 
per unit volume is much larger than in the presently used 
thermal type reactors [10]. The high power density results in 
steep temperature gradients, large coolant temperature rises 
and therefore coolant distribution problems. 

The second important fact is the lack of inherent load 
following capability of this type of plant [X]. This is due 
to the fact that the load fluctuations do not feedback as 
efficiently through the coolant medium as they do in thermal 
type reactors. For these reasons constant automatic super- 
melon dicmcentrol Of a fast breeder reactor is necessary. 

To design a controller for a nuclear power plant eee 


necessary first to obtain a mathematical model. In this report, 


| 
} 


the dynamic equation of a LMFBR are derived. The next step in | 
the design procedure was to determine the dynamic behavior : 
of the complete plant following selected changes in both the 
controlling inputs (reactivity, valve opening) and in the load. 


Because of the nonlinearity of the system, model simulation 


was required in order to obtain a solution for the dynamic 





S (age 


Ae. 


response. From this simulation conclusions were drawn about 
Mme yeharacteristics which the controller must have and about 
the required control strategy. 

The real-time solution of the dynamic equations is performed 
by the analog computer (COMCOR Ci 5000) in a hybrid linkage 
with the digital computer (XDS 9300). 

The controller design developed in this report uses the 
Grgieal Computer as the controlling element. The economic 
advantage of uSing a digital controller results from the pros- 
pect of operating a general purpose type computer in the time- 
Paideemmecde CO perform Doth the control decisions as well as the 
Menltering data processing and such other related functions 
that are presently being performed by such digital computers. 

A digital computer automatic control program is then 
developed on the basis of the results of the simulated pliant 
response. The design strategy includes in addition to the 
ideas of derivative and proportional control, precomputed 
control trajectories. The digital controller design was 
implemented by applying FORTRAN IV programming to the XDS 9300 
Gigital computer. The effectiveness of the strategy for the 
automatic controller is demonstrated by using various step 
changes in the load demand. 

A special software routine was written to automatically 
scalewene analog computer and proved to be quite useful in 
improving the operating efficiency of the hybrid computer 


Simulation. 
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ti ae DESCRIPTION 


The plant considered in this investigation consists of: 
1) Fast Breeder Nuclear reactor with control rods 
2) Primary sodium coolant loop 
3) Heat exchanger (counterflow) 
4) Secondary sodium coolant loop bypass and bypass 
€emenol valve 
SS sveam generator 
6) Steam loop with turbine and condensor 
Additional components which are found in nuclear reactor 
power plants such as pumps, intermediate heat exchangers, 
Pm@essute tanks, equalizer tamis, storage tanks, rek@ater, 
Superheater, and intermediate turbine stages, were not included 
Since they would not. have served the purpose of this investi- 
Pein. By adding these parts, the basic structure of the 
system would not be changed and would still be valid. 


A flow diagram of the a 1s given below. 
Ni! 


CR ee 
er je [ray eee TB3}SGEN 
BOT TBO2 


Fig. 1. Plant flow diagram. 
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Inlet temperature heat exchanger secondary side 


Average temperature heat exchanger secondary side 


ive campcracuLe Steam Generator primary side 
Outlet temperature steam generator primary side 
Average temperature steam generator primary side 
Average temperature steam generator secondary 


side 
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et. MATBEMATICAL MODEL - PLANT DYNAMICS 


A. REACTOR CORE 
1. Neutron Kinetics 
The neutron kinetic equation with one group approxima- 


tion [1] is given as 


nit) ARC) -E a) a eau Ge? (alaliaee) 
de(t) _ B Z 
ae n(t) erates) CEE a2) 
where n(t) = neutron density 
6K(t) = reactivity 
@ = fraction of delayed neutrons 
& = neutron lifetime 
A = average decay constant 
c(t) = average concentration of delayed neutrons 


The steady state conditions are given by 


k(t) = 0 
p a 
Se = Se 
where es steady state neutron density 
Cc, = average Steady setate Coneentration of 


delayed neutrons 








When setting 


B ° 
= reactivity in dollar units 

mit) _ 

— N(t) 

oO 
Gi). 
sean 

O 


the following equations are obtained 





Ge) = 
eee = (N(E) = DUE) AT 


The numerical values for a fast reactor for 8 and 2 are approx- 


imately given as 


B23+ 10° 
eG 10! 
Therefore f 5 - 103 


Hi ehe ditterential equation for N(t) is brought into the form 


dN (t) 
Gite 





Qeeateence) = N(E) © D(t)— = 


The last term, being small compared to the other three terms 
may be neglected. Using this approximation the equation for 


the neutron kinetics in the reactor core are then given as 





GD) _ 

dt en ce) DCE) (TEheo) 
pe, See) 

Nie = ToR(E) (TI 46) 


14 








with the steady state, initial conditions 
N(O) = D(O) = 1 (Claires) 


2. Reactivity Coefficients 
In addition to the controlled reactivity input applied 
MeeeOoSie1tOning Of the control rods, a nuclear reactor has 
inherent reactivity feedback, produced by thermohydraulic 


effects. The total reactivity is given as 


P(e) 


- Ky (t) 1 Ke (ot) (Ite) 


where Kp (t) controlled reactivity input 


inherent reactivity feedback coefficients 


K,, (t) 


In a fast reactor the following reactivity coefficients are 
of importance 

a) Doppler coefficient 

b) Fuel axial expansion coefficient 

c) Fuel radial expansion coefficient 

qd) Clad temperature coefficient 

e) Sodium temperature coefficient 


a block diagram for the stability model is given in Fig. 2. 


1s 
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Sodium temp. 
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Radial Exp. 
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Fig. 2. Block Diagram for reactivity feed back paths. 


iiemmonestcal values of the reactivity coefficients are given 


ace (2 ] 
a) Doppler coefficient ele cents/ % power 
b) Fuel axial expansions ee cents/ % power 
coefficient 
¢) Clad temperature | + .0138 cents/ 3% power 
coefficient 
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a) Sodium temperature + .0057 cents/ % power 
coefficient 
e) Radial expansion - .0043 cents/ 3 power 
coefficient 
From the above data it is seen that coefficients c, d, and e 
are several orders of magnitude smaller than a and b. There- 
fore, these reactivity coefficients will be neglected and it 
will be shown in section IV that this does not lead to any 
noticeable inaccuracy. The reactivity feedback coefficient 


may now be expressed as 


Ka = Ky aE Keyp Chi) 
where Ky = demepler coefficient 
Kop = tuel axial exvansion coetticient 


and the equation for the total reactivity is given as 


Peeeine. th + K 


. R D (Pit) 


EXP 


If the numerical values for the reactivity coefficients are 
introduced along with the dollar units of reactivity (see 
Appendix A), the final equation for the reactivity feedback 


os Given as 


Ko =K, =. Ones. Er Cite ap) 


(een the control block diagram given in Fig. 3. 


a7 
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Fig. 3. Reactivity Feedback Model 
(Samp ila fred), 


3. Heat Generation 
The rate of heat produced per unit volume of fuel is 


eevlead the volumetric thermal source strength Qn and 1S given 


by [4] 
Qn = Gs nt) * of * $ (Tey tye) 
where G = energy per fission, Mev 
n(t) = number of fissionable fuel nuclei/em” 
a - MmerlOsicobi CaliSslon) CELOoSS SECELoOn OL ENG eEnel 
® .= Helieron Llux sper em? Sec. 


ieneic taken aS a constant for a certain fuel element, it is 
seen from the above relation, that the generated heat is 


Preroreional to the product of the fission cross sectagn 


It, 
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and the neutron flux 9. The neutron flux is dependent upon 

the fissionable species used and the energy of fisSion neutrons, 
which in turn, is dependent upon the reactor temperature. If 
Of 1£f determined for a certain reactor, it is seen that the 
heat generated at any point is directly proportional to the 
meutron flux at that point, which is directly proportional to 
the earlier used N(t). The generated heat may therefore be 


written as 
0. ee N (t) ial Tiepiuesy) 


where A = proportionality factor, which is constant if 


Of 1s assumed to be a constant 


Bee abAr TRANSFER 
i. Reactor 
aaa fuel to Co@lmime 

The core of a fast nuclear reactor consists of a 
great number of fuel elements. In a proposed 1 MW Fast 
Breeder Reactor [2] for example, there are 470 fueld rods, each 
fuel section (core length) is 2 feet long and the rod has an 
overall eer Of 25 Incnes. | the wtuc las) enclosed. ned 
@lodeaing of wall thickness of .015 inches. 

Due to the size of the fuel element the radial 
flux distribution within the core may be neglected. The axial 


filme distribution is shown in Figure 4. 








Coolant 


Fig. 4. Reactore Core and Axial 
Flux Distribution. 


iomene axial variation in flux is taken into account, it is 
seen from figure 5 that the generated heat and therefore the 
temperature in the fuel element, cladding and the coolant 


wares in the axial direction of the fuel element. 


U 





Pare 
Elen im: 
Fig. 5. Temperature variation along Fuel 
element (center, cladding, 
COO aime) 
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In this investigation, however, since the reactor 
is only a part of the total power plant, these speclal effects 
will be neglected. Only average temperatures will be 


considered. 





iG. 
O 
_— 







eo eee ater aes er peep 
rene en 
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Fig. 6. Temperaeure profile an 
Eve] element. 


Figure 6 shows a typical temperature profile within 
the fuel element, cladding and coolant together with the average 
temperature profile considered here. The effect of the clad- 
ding is neglected, which is justified in view of the earlier 
given numerical values. Taking the average temperatures 


Simplifies the mathematical solution without introducing 
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considerable errors. If more accuracy is desired, the obtained 
temperature values could always be corrected with a multipli- 
@aeten factor for point values or for optimal accuracy the 
temperature distribution model may Le obtained by the solution 
of partial differential equations [3]. | | 

Under the above assumptions the heat balance for 
a Single fuel element, whose diameter is small with respect to 
Mmiencore, may be wratten for the steady state heat transfer 


eS 6 


Q Ci le 


Qn 2Tr dr L = Ortaars moe 





Fig. 7. Cross section of cylindrial 
fuel element. 
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ft SOlucion to this equation with the boundary condition 





aT _ 
ale, 0 
Tl 


1s (see Appendix B) 


I! 


On KA ULE eee ay) (Tisai 


where K Thermal conductivity, in this case an average 


value for the fuel cladding and coolant 


cal 
2 oO 
sec cm 


A = surface area of the fuel element Lom? ] 
The thermal conductivity of the fuel and ciadding and the 
physical properties of the coolant (density, viscosity, and 
Specific heat) are considered to be constant and independent 
of temperature, which 1S a valid assumption for medium power 
variations of the reactor. For the above assumptions the 
product K . A may be considered constant, Aeseane only upon 
the physical size of the fuel element. For the dynamic case 
of the heat transfer the heat capacitance of the fuel must be 
taken into account. Using the analogy of an electric circuit 
for heat flow, where heat flow is equivalent to current, 
temperature equivalent to voltage, and heat sink (dissipation) 
paphivAlent to a resistor, the following relationship can be 
written | 


eee cue he eter lo) 
N Ee F 
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Fig. 8. Heat flow model. 


The definition of the heat capacitance is [7] 


A QO cal 


gee pees 


iL sec 


Therefore the above equation may be written as 








arr _. i _ ; 
ane = 7 (TAV TF) + A' N 
i KA cal 2 Le il: 
where ee i oa eines aaa 
16 sec cm °C cal 


= fuel time constant 
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ig 
ser aePOr cola lity, £actor 
b. Heat Removal 
The amount of heat removed by the coolant is given 


as 
Of: Ame Ve Cc © AT (EME IETS) 


where 8 = density of coolant Fa 


5 
cm 
A. = cross sectional area of coolant passage [em] 
V = coolant speed| —= 
sec 


C = specific heat of coolant eons 
Gi C 
AT = coolant temperature rise within the element [°C} 
AS shown earlier the temperature rise of the coolant is not 


@ linear function of the position of a unit element within 


Ememcoolant channel. A second nonlinearity accours in the 





coolant channel which is shown in Fag. 92) The temoerature won 
ee 
eae se eee tan 
| : | 
a : WAL | = 
<> | 7 Coolant | TA ae 
CENDRE 
es | > 
| 
| iY 
Cae 7 + YY F LVI S® 
Fuel 


Fig. 9. Temperature in coolant path. 





coolant at the walls of the fuel element .are higher than at the 
center of the coolant channel. This nonlinearity may be 
Mmeglected for the reasons stated earlier, anGd dimaverage! dreetn 
perature profile will be assumed as shown in Fig. 9. The 
specific heat of the coolant and the coolant density will be 
assumed constant. The heat balance equation for the dynamic 
case 1S again derived with the aid of an electrical equivalent 


Beewomown in Fig. 10. 


Coolant 
Re ee NO 


Fuel 





Eigwn O- Electrical equivalent for heat removal 
by coolant from core. 


The heat flow equation is 


_., aTAV ae 
2p = Co “ae * Sav aoe 
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where Qay = Q (EvonwequatLonmiLi, 18) 
which may also be written as 


Qay =mc AT 


mec (TCL - THI) 
coolant flow rate ae 


II 


where m 
sec 


After rearranging and using the relation 


JCA oe Se 


TAV = 5 


to eliminate TH1, the following equation is obtained 


aTAV 1 Ai il i 


























Z 2 2 Z 
where + = KA |_cal com? socaney i] 
eo Co a cna € 7 Se6| 
M—renc gm °C cab. af 
T, ! C, sec Cad come 1 SG 


ae Heat Exchanger 


Assuming no heat loss to the outisde the following 
relation holds: The total change in heat of the secondary 


sodium within the heat exchanger. 


AQ, = AQ, (ele 2am) 


m c(TBI1 - TBO1) 


where AQ, 


the heat transferred from the primary sodium 


Relea TB) 


It 


AQ. 


It 


the heat transferred to the secondary sod™lun 


2a 





Using average temperature of primary and 


secondary Side. 


Equation III.21 becomes now 


Marrs] - TB2) = 


or, after using the relation 





ey = SEeO1 
Z 
TBI1-TB2 = + (TBI1-TB1) 
M1 
Pe 2 hc gm cal sec 
where Ml] =. aa | nae 
If now, as done before, the heat flow 


of the primary sodium and tne overail 


assumed constant, Ml may be seen as a 


m c(TBI1 - TBOl]) 


(Il pee) 


Zac 
en) c -> Jinonet gone 


cal 2 


cm 
rate, heat capacitance 
transfer coefficient is 


heat exchanger constant. 


The dynamic equation for the secondary Side is derived 


as in section II.B.1l.b. 


_ aTB2 
as oe ae es 


where Qp = WA tT) —=eiB2) 
Qo = me (TH2 - TC2) 
Co = heat capacitance of secondary sodium 


Smeaeeer rearranging 





dtTB2_ 4 LL po + & 1 + 
dt T 3 T. T3 
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3 


(Tite) 





where Ll = KA Nec 
T Co 


ub s ZS me 
T3! Cc sec 


Bae Steam Generator 





Since the flow of secondary sodium through the steam 
generator can be changed by the bypass valve in the secondary 
circuit the flow rate of sodium through the steam generator 
will be defined as xm, where x is the valve opening 
Woe~ XxX < (1). 

x = 1 valve fully open (100% flow rate) 

x = 0 valve completely closed 


As in section II.B.2. the following equation is obtained 


iss -— ws) = x m c(TBR2 - TBO2) ({[71I.24) 


eliminating TBO2 and rearranging gives 


ee a 
TB3 - TS = MD (TBI2 TB3) 
Se? eS 
where M2 _ in 


a dimensionless steam generator constant 
under the assumptions of section Tl .BaZ: 
For the dynamic equation of the steamside, the principle used 
in section II.B.1.b. is again applied and the following 
equation obtained 


| os (imneeon 
Q, = Cg ge t Lg | 
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KA (TB3 - TS) 


< 
5 
o 
ay 
© 
1 
it 


Q 
N 


S neat Capacitance Of the steam 


Qq = heat consumed by the steam plant 
From the Carnot cycle it is known, that the work done in a 
thermodynamic machine is proportional to the difference of 


the added and rejected heat. 
Wa (Q, - Qp) 
but Qa Oe Qe 


and Qn is proportional to the entrance temperature to the 


steam plant and Q. is proportional to the steam outlet tem- 


R 


perature. Taking the average temperature it is seen that 


Oe ees 
Rs 


$ 


from this relation Q may be written as 


Q,=K-* y+ TS 
hMe= Constant, assuming, constant flowrate, Nneade 
capacitance and linear temperature chance with. 
reapect. tO PLressure over oe range of operation 
gi 
oC 
Ww = dimensionless proportionality factor Eom ENC 
igadeeaemand (OQ <p < 1) 
) = 1 + 100% load | 


w= 0-7 0 load 


Mereieton (LLlL.25) is then 


= ats 
fuaess) — ES) = Cc. =. + K* y “4s 
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pueeancer rearranging 








— = -(—= + Fo) ee oe + = SS Grr. 26, 
5 5 
meee TX Cal = Cc 
where K' = ee ee dimensionless 
1 = BA TL 
iz Co sec 


er HEAT TRANSPORT 

To obtain dynamic equations which include the time delays 
present within the heat transport system of the plant, the 
Taylor series .expansion will be used. Throughout this section 
it will be assumed, that there is no heat loss during the heat 
transport from one plant component to another. 

ieee nmiary Circuit 

The transport delay, i.e., the time it takes for one 

unit volume element of sodium to get from the exit of the 
meeecror (wath temperature TH1) to the entrance of the heat 
exchanger (still with temperature TH1, but now called TBI1) is 
designated T4. Thus 


apse 


ay at 


fre, LH) (t) CEE. 2%) 


2 a Ge sy 


Expanding the right side of equation (III.27) in a Taylor 


series 
TBI1(t+T,) = TBI1(t) + Ty ompre st) pees (FET, 2) 
where H.O.T. = higher order terms. 


3h 





Assuming slow temperature changes, which is reasonable 
Bor a large stationary power plant, all higher order terms may 
be neglected. 

Replacing TH1 and rearranging, equation (III.27) may 


be written as 


ait Ly Ty v4 





For the transport delay from the heat exchanger back to the 
reactor the same time duration will be assumed. 
Therefore TBOl = TCl 


me SBOU(t) 


EGC Ty) CEI .38) 


Using again Taylor series expansion, replacing TBO1 and 
rearranging, the dynamic equation for the heat transport from 


the heat exchanger back to the reactor is given as 


ane 1 2 1 
== = - = TC] + = TBl - = 
dt Ty, a i 


Pet Gt tee; 1) 
ee vcecondary Loop 


The heat flow diagram for the Secondary sodium loop is 


~ shown in Figure 11 where Coe (N =si 2 DAS Pile. alec mene 
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Pig. ll. Heat flow diagram for secondary 
Giccuit, 


which heat is transported by a unit volume element. The 
function of the bypass valve X 1s to increase or decrease the 
amount of sodium flow through the steam generator: it does 


met influence the temperatures. This may be shown as follows: 


usince the equivalent of an electric circuit, the "model 


equation" for node A is written: 


a ts (ELE ee) 
eee 3 


aneec TBI2 =m c TH2 - (lL —x)m ¢ Vie 


SLB LZ 


eZ a ee) EZ 


MB Z 


DHZ 


Sis 





Bem orore the same procedure aS in section III.C.1, may be 


used to arrive at a dynamic equation for the heat transport 


G@etay from the heat exchanger to the steam generator. The 


Pimemcauration of the transport delay wall be called T_. 


Since TH2 PE 


and TH2 (t) Peay Ce +L 


5) 


and, again using Taylor series expansion for TBI2(t + Te) and 


replacing TH2 (t) gives 


ATBI2(t) _ 
dt 


oi 


TBI2 + < TB2 eee (IIL. 


5)5)), 


For the heat transport from steam generator to heat exchange 


the "“nodel equation" for "node" B is 


on 0, + Q) Cum 


me TCc2 = xm ce TBO2 + (1l-x) m ce TH2 


Assuming that the duration of the transport delay is again 
for the flow through the bypass and for the flow from the 


steam generator, equation (III.34) becomes 


POZE + T.) = x TBO2Z(t) + (1-x) TH2(t) (ie 


Replacing TBO2 and TH2, using the relations 


M2 _ TBI2 - TB3 
x iM — TS 
TB3 = Tei2 + TBO2 
Z 

Z 
TB2 = he * Ed 
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eaemexpanding TC2(t + Te) in a Taylor series, equation (III.35) 





becomes 
arc2_ oe wa i 
See oto oe * mon 8B2 + aa TS + aay PBI2 (III. 36) 
5 5 5 5 
where 
p os aie 
5 2-X 
PP 25 
5 2-2X 
T_ (x+M2). 
7am sD 
5 2x M2 
T, (xt+M2) 
ae Se 
5 x (X-M2) 
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IV. SYSTEM SIMULATION 


A. ANALOG COMPUTER PROGRAMMING 

The coefficients of the differential equations developed 
mieoeceton If were calculated for a steady state solution (see 
Appendix C). Using the numerical values of these coefficients 
the following set of equations which describe the dynamic 


behaviour of the plant was obtained: 


D = £.58(N-D) 
TF = - .5 TF + .5 TAV + 725 N 
TAV = - .5333 TAV + .5 TF + 4.8333 TCl 
TCL = - .25 TCl - .0988 TBI1 + .3488 TB2 
PBIi = - .25 TBIl + .5 TAV - .25 TCl 
TB2 = -6.9626 TB2 + 3.0233 TBI1 + 3.9394 TC2 
Tene = - .25°9RI2 +S TRQee 25 TCA 
me2 = —- ,2625 TC2 + .275 TS - .0375 TB2 + .025 TB2 
T = -—- .0421 TS - .0168 + p + TS + .0421 TBI2 
N a 


D / (1-K,) | 
To simulate these equations on an analog computer and to obtain 
Naeem accuracy, the gain coefficients eee to be scaled with 
respect to the maximum voltage level of the analog computer 
(see Appendix D). The scaled equation used for the simulation 
are given below. 

| D =  .58 (N-D) 


- ,5 TF + .1833 TAV + .4833 N 


TP 


=5e3230 TAV + 1.3636 TV + 4.3937 30er 


TAV 
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I 
t 


mel Pople = 21235 BBli +. 23214 @e2 


TBIl = - .25 TBI] + .44 TAV ~ .2 TCl 
TC2 = — .2625 TC2 + .2895 TS - .0474 TBI2 + .025 TB2 
TB2 = -6.9626 TB2 + 3. + 78 TBI] + 3.9394 TC2 
TBI2 = - .25 TBI2 + .3958 TB2 - .1979 TC2 
TS = - .0421 TS - .0168 « 4 + TS + .0505 TBI2 
N = D/(1-K,) 


The temperature variables for example are no longer in °F but 
in meres . 

A block diagram, using both the symbols shown below and a 
linear flowgraph for representing the analog computer setup 
may now be drawn (Figure 13, 14). (See Appendix E for the 
analog computer diagram corresponding to the flowgraphs shown 


in Figure 14) 


+> Adder Amplifier 
4\>— Integrator 


{)- Potentiometer 
4+<|- Multiplier 


ot D1 Vider 


Figure 12. Symbols for analog computer components 


It proved very valuable to write a simple Fortran executive 


program which performed the above calculations each time a 


a7 








Fig. 13. Analog Computer setup -— reactor core. 
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test or simulation run was made. This program contained the 


following steps: 


a) 


b) 


ec) 


d) 


e) 


=) 


g) 


h) 


Hie) 


Read in steady state plant variable values, 
constants and chosen maximum plant variable 
values 

Calculate the numerical value for the coefficients 
of the differential equations 

Calculate magnitude scaling factors 

Calculate scaled coefficients for the analog 
computer components 

Gateuracte initial conditions for analog 
integrators 

Write out all calculated values and the scaled 
and unscaled differential equations 

Calculate the initial steady-state output vaiues 
of the amplifiers 

Set all potentiometers 

Scan all amplifiers and compare with calculated 


values, write out errors 


The amount of time saved by uSing this executive program 


during the many runs performed was appreciable. The program 


had the following advantages: 


a) 


b) 


Equipment failures were detected immediately and 


isolated to only a few components by the test part of 


the program. 


The scaling procedure was done by the computer. Since 


the maximum plant variable values were not known, the 


optimum scaling coefficients had to be found by trial 
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and error in many runs. By just changing one statement 
in the program (maximum variable value) the analog 
computer was ready for the next run within seconds. 

c) The documentation was done automatically by the digital 


Computer. 


Pe DYNAMIC CHECK 

A theoretical curve describing the behaviour of the reactor 
of the proposed plant following an accidental reactivity input 
was available from Ref. [2]. The validity of the dynamic 
Sweatious describing the reactor was demonstrated with the aid 
@eetnis curve. 

To check oe solution of the analog computer the system 
equations were also solved on a digital computer (IBM 360/67) 
memng 2a Lourth order Runge Kutta routine with accuracy testing 
eeecingle and double increment and double precision calculation. 
The results of these. dynamic checks are stated below. 

ime comerolled Accidental Reactivity Input 

The accident considered assumes that two control rods 
Grop into the core. This would correspond to a reactivity 
insertion of 100 S$/sec. It is further assumed that scram 
Action takes place .2 sec. after the accident starts. 

A solution of the system equation to this accident was 
obtained on the digital computer with the above described 


mencine. 





Because of the nonlinearities in the plant dynamic equa- 
euons, an analytical solution for the calculation of tne 
maximum values cannot readily be obtained. 


Al 





The digital computer solution gave exactly the same 
results for the temperature variables TAV and TF. The numerical 
values for the neutron density N(t) varied slightly in magni- 
tude and in the time of occurance. This slight difference can 
be explained by the simplifying assumptions which were made in 
deriving equation I1I.5 and II.6. The one group delayed 
neutron and the prompt jump approximation effects the time of 
occurance of the maximum value of N(t). By neglecting small 
reactivity coefficient (equation II.10) the magnitude of N(t) 
is influenced. In spite of this small difference in the 
result, the equations used for the simulation of the reactor 
core are thought to be accurate, especially since, as stated 
above, the temperature TAV and TF followed very precisely the 
given data. 

2. Change in Load 

The first part of the check on @he aceuracy of the 
analog computer solution was made by obtaining a digital compu- 
ter solution to a load change in the form of a step change of 
magnitude - 20% for 20 sec. The results are given in Table I 
and in Figure 15. 

| See Chonda in Reactivity 

A second digital computer check of the analog computer 
wimtton consisted in introducing a disturbance in the reac- 
mevety input in the form of a +20 percent reactivity step 
change for 20 sec. The results are given in Table II and in 


maeoure 16. 
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Tiame(sec) TCI(°F) TB2(°F) Totes) 


0 800-2 885.1 750.2 
2 SOUlaeeCsG.6 686 762-3 
4 803.9 890.4 773.1 
6 808.2 895.5 782.9 
g Ces acos 791-6 
10 820.2 908.6 800.2 
15 $35.3 924.9 9818.5 
20 840.1 941.5 834.7 
22 844.6 946.5 827.6 
24 $52.2 948.5 9821.7 
26 857.4 948.8 816.9 
28 858.3 948.1 812.9 
30 860.6 946.8 809.2 
35 855.2 942.2 804.0 
40 850.4 937.3 800.0 
TABLE I 


Digital computer results; temperature changes due to 
Seep Ciange in load of -20% for 0 < t < 20 sec. 
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Fig. 15. Analog computer results Load change as in table I. 
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Taimetsec) TECE(°F) TB2 (°F) SCE} 


0 800.1 OG. OZ 
2 S00. 3 Sod, 6 Tou. 2 
4 Sin’. 1 eo kaG 1 30..5 
6 eZ. 2 CoS vei. 0 
8 803.3 GU 5e5 vole. 8 
10 804.3 oor. Di 260 
aS $06.5 899.64 153.6 
20 808.9 201.8 woe .9 
22 BOIe/ 0G. © 760.2 
24 S070 Coen o VoL. 0 
26 | 810.0 896.4 Ger) 
28 810.0 SUSnI6 162-0 
30 810.0 S923 (IDA @ 
35 810.4 Oa. 2 Wo 4 
40 S206 Bp. 0 Four > 
TABLE II 


Digital computer results; temperature changes due 
to step change in reactivity input of + 20 cents 
moro =< tt < 20 sec. 
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Fig. 16. Analog computer results Reactivity change as 
in table II. | 
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4, Discussion of the Dynamic Check Results 


It can be observed from the above results that the 
developed equations are valid for the plant under consideration. 
Comparison of the given digital computer solution and the 
analog computer graphs shows that the analog computer gives 4 
solution which is within the bound of less then one percent 
exact (+ .5%). 

The results of the tests have demonstrated the validity 
of the simulation within the above bound of accuracy and under 


the earlier made assumptions. 


47 


VI. CONTROLLER 


A. CONTROL STRATEGY 
The inherent incapability of the plant to follow load 
changes which is due to the small eeaeenae feedback coef- 
ficient of the coolant, is a disadvantage in the control of a 
Fast Breeder Reactor. It makes complete control of the coolant 
temperature in the reactor necessary. An automatic controller 
for power level control must be developed whose two main 
purposes are 
a) coolant temperature level must be kept constant 
to obtain optimum breeding ratio and optimum perfor- 
namee of “the reactor and 
b) steam conditions are to be kept constant for optimum 
performance of the steam plant. 
On the other hand several advantages which are inherent 
for the plant help the controlling action to achieve smooth 
and efficient power level changes. The reactor is. automati- 
cally isolated from the load by the secondary sodium loop and 
to a certain degree by the primary sodium loop. Therefore, 
the inlet temperature to the core changes only very slowly. 
These two circuits also serve as a heat reservoir. This allows 
a temporarily larger rate of load change than the rate of power 
Change which the reactor can achieve. 
The rate of power change obtained from the reactor is 


memaeced by the maximum rate of change of reactivity. This 
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. cents 
value is taken as 2 


J oeer ile |) heen clconlamieed by 
mechanical factors such as thermal stress within the reactor, 
heat exchanger or ete generator and speed of mechanically 
driven valves. The thermal stresses were not investigated 
here, it can, however, be stated that in none of the test runs 
mee limiting values for the rate of change of temperature, 
e.g. the coolant in the core with 100°F/sec., as given in 
Ref. [2], was observed. The limiting factor in this investi- 
gation was found to be the rate of change of reactivity. 

The control strategy is based on the following thoughts: 

Because of the extreme long time constant, a pure pro- 
portional controller, which acts only following the occurance 
fm EtOns 1m the to-be-Ccontrolled quantities was not sufficient. 
Therefore, an anticipatory or derivative control had to be 
mimebemented. To get the best derivative control action signal, 
meewderivative or the rate of change of Corpor atures which had 
to lead to changes in the to-be-controlled temperature had to 
be sensed at the earliest possible time. This steady state - 
memeat most very slow power level change - strategy is aoe 
valid for larger power level changes. In this case preprogrammed 


fenerolling action had to be implemented. 


Se CONTROLLER DESIGN 

The choice was given between a pure analog-logic controller 
@e a digital controlling program. 

The analog controller had the advantage that it could have 
been implemented comparitively easy for a plant simulation of 


meme cize, On the other hand it was felt to be more realistic 
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to design an automatic controller which uses the digital 
computer since modern power plants are equipped with digital 
computers used for many purposes - from bookkeeping of the 
plant datas to paycheck writing - that a digital controller 
would represent only a small and comparitively low cost 
addition to the main computer. Investigations have also shown 
for, “Mat for large scale plant controllers digital controllers 
Become more economical than analog controllers. Therefore, it 
was decided to design a digital controller. The following 
controlling devices were "built" in to the plant: 

a) control rods. It was assumed that the rate of chance 
of reactivity could be varied. For the preprogrammed 
controlling action the maximum rate of 2 cents/sec. was 
assumed; for the steady state controlling action a 
value of .1 times the maximum rate was found by 
experiment to give the smoothest control dynamics. 

by Bypass valve in secondary circuit. The opening time 
(with maximum speed) for the preprogrammed phase was 
assumed to be 20 sec. from completely closed to fully 
epem, fer the small controlwane action .2 of this walwe 
was used. 

Because of the constantly changing value of the bypass valve 
opening, several gain factors in the dynamic equations 
developed in section II changed during operation. As a result 
of this the analog computer setup had to be modified (see 


Appendix E). 
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The above described controlling devices where controlled 


by the digital computer which in turn received the data from 


the simulated plant on the 





Control 


nage? selant and Control! points. 


Emalog computer (see Fig. 17). 
Eines Conerol strategy called for 
Eveee@ee programmed power level change control. The control 
trajectory for this control was obtained by ie GG ein 
ieacmemanges and correcting the errors in TCl and TS 


Cl Ne ' 


Wien manually operated potentiometer, which simulated 


Syl 





control rods and bypass valve opening. The resulting 
datas were evaluated with the aid of aleast square 
straight line approximation program to give the 

Gemerol tragjeccory for the” control reactivity and valve 
opening as a function of the load (see Appendix F). 

b) Derivative control. 

Following the thought described in IV.1, the rate of 
change in TAV and TB2 were calculated by storing values 
every Mou of a second and then calculating the rate of 
Glipmge over One Second. 

e) Proportional control. 

The to-be-controlled temperatures TCl and TS were 
compared every ines of a second with the set value 
and the error was calculated. 

To take the addition of noise in actual plant measured 
variables into consideration, all measured values were put into 
a dead zone. This had also the advantage of smoothing the 
@entrolling action. 

A flow chart of the basic program for the digital controller 
ieeglven in Figure 18. The individual parts of the program 


are described in detail in Appendix G. 


eZ 





Start 


Scan Load 


4 
| f 


Has 
No Load 
Changed 
Yes 
Pre programmed | aoe | 
Controller "| wernt yo 
eactivity = F(Load) aa | 
Bypass = F(Load) oo 
achieved 


Scanload 





Small Deviation 
Controller 
Proportional - 
Derivative Control 


Pig tO. whlow Chartster digital controller. 
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VI. RESULTS AND COMMENTS 


The results which were achieved by controlling the plant 
Meee the described controller are shown in this section as 
stripchart recordings in Figure 19 through Figure 21. The 
test inputs were step changes in load of different sizes and 
pe GQiftferent time intervals. 

Figure 19 shows the plant dynamics following a step change 
in load of -50%. This step change does not seem to be very 
realistic, it proves, however, that the initial overshoot in 
the to-be-controlled temperatures of little more than 100°F 
is damped down very effectively to a secondary overshoot - or 
first undershoot - of about 20°F. The initial overshoot 
cannot be avoided since the controlling devices work at 
Seria allowable speed during the preprogrammed part. The 
Same statements hold for the step change in load of +50% which 
was introduced after the transients had almost settled. Even 
a load change of this amount does not endanger Eee Paanee Une 
to the fast working and direct preprogrammed controller - which 
may be seen from the value of N which increases only saboucsiZ- 

It is believed, that this test input and the resulting 
Meanszent Gemonstrate the validity and necessity of the 
feyeslLon of the control program into its two main parts, it 
m@lso demonstrates effectiveness of the preprogrammed controller. 

Figure 20 shows the plant dynamics following a step 


change of different magnitudes and at time intervals Which 
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are shorter than the settling time of the dynamics. A 20% 

Step change in load produces only a 25°F, or 35°F respectively, 
Syvershoot in the controlled temperature. The increase in N 

is only 4%. The interplay of prepregrammed and small deviation 
control may also be seen from Figure 20 by considering the 
reactivity and bypass valve curve. 

Figure 21 was taken over a time of more than 1 hour. 

Though step changes in load of 10% do not seem to be very 
practical this recording shows a more realistic picture of the 
plant controller. It especially proves the accuracy of the 
Steady state controller. After the settling of the dynamics 
this controller is able to keep the error in the controlled 
menperatures within a bound of less than + 3°F at a controlled 
temperature level of 800°F for TCI and 750°F for TS. 

The above results demonstrate the effectiveness of the 
designed controller for the plant layout simulated in this 
investigation. Since conventional stability calculations like 
meot-locus plot, Bode plot or Nyquist diagram cannot be 
applied directly to a controller of the kind deve loped nere, 
this topic is not included in this report. It can however be 
Stated that no instabilities following any of the applied load 
changes were observed during the simulation. 

The addition of more plant components in the Simulation 
of the system will, as stated in section IIT, not change the 
basic form of the dynamic system equation. It will however, 
change the dynamic response to a certain degree. The same 
is true for the nonlinearities which were not taken into account 


in the simulation in this report. 


22 


Further investigations can include the above considera- 
mpons, however, a simulation of a complete plant which includes 
all components and takes all nonlinearities into account is 
Only possible on a very large scale digital, or better, hybrid 
computer. If a computer of this size is not available an 
accurate Simulation can be done component by component as in 
Ref. [3] for a LMFBR. Accurate stability investigation may 


then be performed by simulation. 
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APPENDIX A 


m CALCULATION OF THE REACTIVITY FEEDBACK COEFFICIENTS 
The reactivity coefficients considered are 
a) Doppler coefficient 


b) Fuel expansion coefficient 


a) The Doppler coefficient is a nonlinear function 


of the fuel temperature and is shown in Fig. 22. 





oO 
Be 


1000 000 000 4000 


Pidamcc. Moppler Coctircient, 


Since the operating temperature of the reactor is at 2400°F an 
average doppler coefficient of SAeceaon could be taken. It 


can be seen, that even at temperature variations of +400°, 
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mach have to be avoided in a fast reactor, the change in 


@eppler coefficient is at most ee Omar TO accom: for a 


certain degrees of uncertainty a conservative numerical value 


.3.10°> 


will be used. Changing this value to dollar units 
of reactivity with 6 taken as 29.88x10"* the doppler coefficient 
1s obtained as 


on 5 | 
re Sees eo” ~=S aS Lars 


a =A 
2976610 


b) The fuel expansion coefficient is considerably 
smaller than the doppler coefficient and since its reliability 


for axial fuels is highly uncertain the conservative numerical 


6 


value of -1.0°10 ° will be taken. The conversion into dollar 


Paes reactivity gives: 


ie ome 


4 


K = 2. 834500 yeas 


EXP i‘ 
29 .88x10 


The total feedback reactivity is therefore given as 


Ke = Rp tr Koy 


=20004> - .<00034 


eS 6 


-.0014 dollars 
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APPENDIX B 


me oOLUUTTLON OF THE STEADY STAIL, SECOND ORDER DIFFERENTIAL 
EQUATION FOR THE RADIAL HEAT FLOW IN A CYLINDRICAL FUEL 
ELEMENT 


The heat balance equation is given as 


QO CB aly) 


Q 27x a — Ore FeO 


The differential heat quantities 


dQ 
_ o ig 
ordre Q. a oo — Oy T ay dr Len) 
ian Ge ar 
Q.. = Ky A ae Ke Zip ie le ne (Bes) 
equation (B.2) becomes now 
aT a°T ar 
ota = 277 L KL. aR = 77) Ge Khe FZ i ana. Glia (B. 4) 
ar 
or 
aca 
QO, 27Tr dr L= - 27 L K(r —=+ =) ar (BiS)) 
F ar? dr 


which, after concellation and rearranging, may be written as 


Z 
aT K 
a>” Yr 


aT be 
: ee Qn = ( (BG) 
r 


This differential equation can be solved by standard methods, 


the solution is 


Besa ep Coen xr F.C 


mes Fr aK (1 


> (B.7) 
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The boundary conditions are 


aT 


I 
oO 





ee) 
il 
bi 
KF 
Il 
FY 


r=0 r=R 





QR 
See > aK 
Since 
ae LO (B.8) 
R iE ' 
me f£OlLiow that 
a= TR- L Q (B.9) 
S E : 


Pepeecomorning (B.8) and (B.9) 


Q, = 4m K L(T, - Tg) (B.10) 
Setting 27 RL=A, the surface of the fuel element 
To 7 ts 
oe > Z kK. A —_— (ey) 


Equation (B.11) holds for the region inside the fuel. For 
the cladding region KeaaaG would have to be used instead of Ke 


and for the region in the coolant K Since am@rage 


Coo lant. 
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temperature are used and in addition the individual thermal 
conductivities are considered constant, an average thermal 


conductivity will be taken and (B.3) may be written as 


a 6k 
oF = Ree Tr qr (B.12) 
where r= the distance measured from the center of the fuel 


element. 
To 1s taken as TF, the average fuel temperature and To 
as TAV, the average coolant temperature. 


Equation (B.5) becomes then 


Q, = KA (TF - TAV) (B.13) 
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APPENDIX C 

ee ee UCUGATION OF THE GAIN COEFFICIENTS OF "THE SYSTEMS 

DIFFERENTIAL EQUATIONS 

The gain coefficients in the individual equations may be 
solved for by using the numerical values given for a specific 
plant. This would require a certain amount of engineering 
peeroximation, Since exact values are not available for all 
quantities. Instead of following this approach, only given 
steady state temperature levels for 100% power (W = 1) ata 
bypass valve opening of 953 (X = .95) and given time constants 
were used. These values are [9]: 


a) Steady state temperatures 


LY = 2400°F 
THL = 1100°F = TBI1l 
Tee = 800°F = TBO] 
Tae) = 1050°F = TBi2 
TE2 = 720°F = TBO2 
aS = fs0-F 

b) Time constants 
i = 2 sec = T2 
BS = 4 sec = T4 = T6 
iS = 10 sec 
7 = .1l sec 


Peqguewonal constants used [3]: 
>» = .58 


8 ~0029 35 


il 
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Setting these constants into the obtained dynamic equation 
these equations were "forced" into an accurate steady state 
solution by solving for the unknown gain coefficient in the 


steady state form of the differential equation (= = 0). 
This procedure is demonstrated with one example: 


The equation for heat transfer within the core was obtained 


ass 


opp = - — op + — pay + A'N 
a “ 


which becomes in its steady state form 


(2400) + = (950) + A’ (2) 


on) 

il 

] 
No] 


and soiving for A‘ gives 


= 70 


which is the desired coefficient for N. Using this approach 
pomeall equations individually leads to the system equations 


given in section II. 
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APPENDIX D 


=  oCALING OF THE SYSTEM EQUATIONS 

As mentioned in section III the actual equations must be 
scaled with respect to the maximum voltage levels of the 
analog computer. It is desirable to get maximum output values 
of the equation variables close to the computer maximum 
voltage level to maintain maximum accuracy in the solution of 
piewetecerential equations. To be able to perform this scaling 
in an optimum way an exact knowledge of the maximum values 
of the plant variables 1s necessary. If these variables are 
not known exactly, optimum scaling is done by trial and error: 
a maximum value is assumed, the equations are scaled for a 
computer run and the output signal 1S observed. Etech s 
output signal is either to high or too low, the equations have 
me be rescaled. 

The magnitude scaling procedure uses scaling coefficients 
which ‘convert given plant variables with their individual | 
Pe Giicns ate computer variables with the deminsion volt. 


Taking the example of TAV it is written 


TAV = SLTAV : TAY 


where SLTAV 


scaling coefficient for the plant variable 


-{[ Q@ 
TAV with dimension ei 





PAV = computer variable for average coolant temper- 


ature in core with dimension (V) 
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itewanalog computer used here (COMCOR Ci 5000) had a 
maximum voltage level of + 100V. Assuming a maximum value of 


1100°F for TAV the scaling coefficient is obtained as 


Using this procedure the following numerical values for the 


individual scaling coefficient proved to be best suited: 


SLD = 02 
SLTF = 30. 
SLTAV = 1l. 
SLTCl = 10. 
STBIL = 12.5 
SLTB2 = 9.5 
SIBI2 = 12. 
Shbc2 = ore 
SLTS = 10. 


The scaling procedure for a single dynamic equation 
proceeds as follows: 


Taking equation (I1.20) 


TAV = 5.3333 TAV + .5 ' TF + 4.8333 TCl 


Replacing the plant variables by computer variables gives 





eee 8 6hTAV = - 5.3333°SLTAV* TAY +... 5*+6b2r + iF 


a 


+ 4G 30ne ole ker 


Dividing through by SLTAV and using numerical values gives 


tev = - 5.3333 TAV + 1.3636 TF + 4.3939 ow 


a 





Peg thas procedure for all equations and steady state 
meameonciwves the set of scaledWequation given in JII and the 
Peeresponding initial conditions for the integrators of the 


analog computer. These initial conditions were found to be 








ae (O) i=. ..80% (V) 
jay (0). = 86.3636 (Vy) 
TCL (0) = 80 (V) 
TBI1 (0) = 88 (V) 
mee (0) = 93.1579 (W) 
mem(0) =. 87.5 (Vv) 
mee (0) = 75.7895 (Vy) 

Te O)ee=. 7 5e (Vv) 


The steady state control rod reactivity required is found 
to be 3.36 dollar which, after being scaled, ccerresponds to 


Permo (Vv). 
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APPENDIX E 


A. ANALOG COMPUTER SETUP FOR CONTROLLER IMPLEMENTATION 

It was shown in II, that several gain factors in the system 
dynamic equation were a function of the bypass valve opening 
mR. Since the bypass is used as one of the controlling ele- 
ments in the plant and therefore continually changes its 
numerical value, a way had to be found to continually change 
the gain factors which are functions of X. For the simulation 
without controller the gain factors were implemented as com- 
parecer Set MeeeKevonccers, which take a considerable amount of 
time to be set or adjusted. A continuous dynamic solution or 
Simulation of the plant was therefore not possibie with tis 
setup. | 

The following method was used to get continuously changing 
gain factors and therefore a continuous Simulation: 

The potentiometers were replaced by electronic multipliers 
which had as input the to be multiplied variable (temperature) 
and the euesut of a trunkline which was controlled by the 


digital computer. Figure 23 


TS ye TS K «w-TS 
a —————————— a 


EOw TRUNKLINE 
peu Le set toK [Vv] 
y : 
+100v 


Fig. 23. Implementation of varying gain 
factors. 
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shows the analog implementation for the second gain factor 
mimes Wttferentialsequation for TS (K**TS). The voltage 
representing TS is multiplied by a voltage representing yp 
momeyi < wv < 1OO0/V) corresponds to 0 < wp < 100% load). wv 
was changed by a manually set potentiometer on the analog 
beard. The output of this multiplier is fed into another 
multiplier with the trunkline as the second input. Each 
time (scan intervall = .1 sec) the controller changed X, it 
Simultaneously calculated all dependent gain factors and 
adjusted the corresponding trunkline. Due to the electronic 
multiplication continuous adjustment of the gain factors with- 
out any Me lag was thereby achieved. | 

Figure 24 shows the resulting analog computer setup for 
the plant part without the reactor. For more clarity of the 
graph no numerical values were set in and all initial condi- 
tion inputs to the integrators with the corresponding 
potentiometers were left out. Figure 24 corresponds to the 
Jinear flowgraph of the plant system equation given in III. 

Since there are no changing gain preteen: in the reactor 


Meme of the simulation Figure 13 ian III issnet changed. 
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Xt-1X 


- Fig. 24. Analog computer setup for plant. 
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APPENDIX F 
eee CONTROL TRAJECTORY FOR PREPROGRAMMED POWER LEVEL 

CHANGE CONTROL 

300 data points were obtained by changing the load over 
eee waole range from 100% to 0% in small steps. For these 
rams the analog computer was in 1000 : 1 timescale mode, 
therefore the steady state error for a certain load was 
immediately available. (see Figure 27). Using hand set 
potentiometers voltages ‘were supplied to 

a) the analog computer which corresponded to control 

Meactruyaty 
b) the digital computer which corresponded to the 
bypass vaive position. 
The digital computer calculated the variable gain factors and 
adjusted the trunlines. 

The "control reactivity voltage" and the "bypass voltage" 
were adjusted until the steady state error was zero. The 
obtained data points were then used to obtain, with the aid 
of a least square polonomial fit program, straight line 
approximations as precomputed control trajectory for load 
changes in the plant. These data points are shown in Figure 
Beemer Figure 26. 

As may be seen from Figure 27, it is, with the controlling 
aids of this plant simulation possible to cemtrol TCI 
@empletely over the whole load range, TS can only be controlled 


over a load range from 100% load to about 20% load. The stecay 
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state error for TS at zero load is about +80°F. This is not 


feeestriction on the controllability of the plant. 
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Fig. 25. Experimental data points. 
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Fig. 26, Experimental Data points. 
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APPENDIX G 


A. DIGITAL CONTROL PROGRAM 


i. 


General 


The scan interval was chosen as .1 sec.,i.e., the 


digital computer went through the designed control loop every 


.1 sec. and performed the necessary Comerol action by “scanming 


the assigned variables, performing the necessary calculations 


to get the derivatives or errors, and then going through all 


the logic steps, implemented as Fortran IF-statements. If 


control action was necessary it was performed by calling the 


following subroutines. 


a) 


REACP: Increase reactivity at the rate of 2 cents/ 
sec. Because of the scan interval the voltage of the 
trunkline which represented reactivity input was 
increased by an amount which corresponded to .2 cents/ 
.l sec. Inspite of being a step change the short 
time interval in which these changes were performed 
made the controlling action an almost continuous one. 
The time delay of an electric motor which would be 
used to perform this “controliing ection in a real 
plant, was neglected. This 1s justifiable in view 

of the short time constant of a motor with respect 

to the time constants of the plant. 


REACH : Decrease reactivity at a rate of 2 cents/sec. 
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APPENDIX G 


A. DIGITAL CONTROL PROGRAM 

1. General 

The scan interval was chosen as .1 sec.,1.e., the 
Gigital computer went through the designed control loop every 
.1 sec. and performed the necessary control action by scanning 
the assigned variables, performing the necessary calculations 
to get the derivatives or errors, and then going through all 
the logic steps, implemented as Fortran iF-statements. ete 
control action was necessary it was performed by calling the 
following subroutines. 

a) REACP: Increase reactivity at the rate of 2 cents/ 
sec. Because of the scan interval the voltage of the 
trunkline which represented reactivity input was 
increased by an amount which corresponded to .2 cents/ 
Peace iis pace OL being a Step change the short 
Piteminccmuan Im whwen tmese changes were performed 
made the controlling action an almost continuous one. 
The time delay of an electric motor which would be 
used to perform this controlling action ina real 
plant, was neglected. This is justifiable in view 
of the short time constant of a motor with respect 
E@OmeneCmEame Cemstants of the plant. 


b) “REACM: Decrease reactivity at a rate of 2 cents/sec. 


8 i 





C) 


d) 


oy 


i) 


The 
a) 
D) 
Cc) 
These 3 


ae 


RACP and RACM: same as under a) and b) but with 


the rate diverted by 10. 


Subroutine REACP and REACM were called if controlling 
action was performed by preprogrammed controller. 
Subroutine RACP and RACM were called if controlling 
action was performed by small deviation controller. 
VENTP : Open bypass valve at described rate (0 - 100% 
in 20 sec.). The subroutine changed the value of X, 
calculated the new gain factor and adjusted the 
corresponding trunklines to the electronic multiplierxs. 
The same statements made in a) with respect to 


continuous action and time constant hold also here. 


VENTM: Decrease bypass valve opening at the described 
awe . 
VETP and VETM: same as under d) and e) but with 


rate divided by 5. 


Subroutine VENTP and VENTM were called by the pre- 
programmed controller subroutine VETP and VETM by the 
small deviation controller. 

digital controller program was divided in 3 parts: 
Routing 

Preprogrammed controller 

Small deviation or steady state controller. 


parts are described below. 


Routing 


The decision whether to go into the preprogrammed 


Paee OL tie controller or to the small deviation controller 
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was based on the rate of change of the load. For a load 
change of less than .2%/sec. the control action was performed 
by the small deviation controller. The flowchart of the 
Coiling part: iseagmvyen in Fig. Doi 

Saye eeprogranmmed controller 

The purpose of this controller was, to bring the 
controlling devices directly to their set points, i.e. to the 
points which were calculated by the straight line approximation 
=eas ae tunction of the load - to lie on the control trajectory. 
This controller performed its function without regard of the 
dynamics of the plant. The load was scanned, the set Dorm 
calculated and the respective subprogram (Controlling action) 
@eomered. Fig. 29 shows a flow chart of this part of the 
pereitea CoOMiine il lex. 

Moma Deviation or Steady State Controller 

As stated in III this part of the controller used 
proportional and ace control. The decisions of the 
digital computer were based on the logic tables given as 
Table III and Table IV. The procedure used here is the same 
wee nN analog, logic circuit, instead of writing expressions in 
Bolean algebra, Fortran IF-statements were written for the 
fegee Gectstons of the digital computer. 

Controlling action for change of reactivity was based 
em oay amd ETC] (error in TCl). The logic table used is 
Table III. 

Controlling action for change of bypass valve value 


was based on TB2 and ETS. The logic table used is Table IV. 
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With respect to earlier thoughts any decision made on 
the basis of Table III was weighted higher and if controlling 
action tas required, Table IV was neglected, i.e. the controller 
Pereccerarst £6 bagiag ETE® to zero and keep it zero before ETS 
was controlled. This was done because of the dynamic of the 
plant. No logic scheme could be arrived at, which would have 
made combined control possible. It will be shown in section 
V, that the performance of this controller was very good. A 
flow diagram of this part of the digital controller is given 


a Ergure 30. 
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Start 


Array Load (I) 


I=1,10 
Scan Freq, 10 cps 








Load (I) = Load (I+1) 


Scan Freq, 10 cps 


Scan Load (10) 


& Load = 





Load (10) - Load (i) | 
asec < A Load 


ae, Sec 


Preprogrammed 
Controller 


Small Deviation 
Controller : 


Fig. 28. Routing part of digital controller. 













Routing Part 


of Controller 





Reactivity 


RvR Ss = F (Load) 


Call subroutine for 
necessary control 





action 
Scan Load | 
\ | 
is 
set Point No 
+ Scam Freq, = 10 cps 
‘Achieved 
y 
Yes 


f 


Fig. 29. Preprogrammed controller part of digital controller, 
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Logic Tables for Small Deviation Controller 


TAV ee: "°° Reactiv. 
0 0 7 
a; 0 D 
_ 0 I 
0 + D 
+ + D 
- </> 
0 a I 
+ - a 
~ - I 


Table IIL. Re@etivity control. 


Eye Eas Boel vaLve 
0 0 sts 

+ 0 D 

- 0 iE 

0 + D 

ats ate D 

= t Pyar 

0 = iD 

+ - w/e 

= os I 


Table IV. Bypass valve control. 


Symbols: 
+ EGEOmw Or Dery. > 
- PrerOnmoi DELrLV. < 
0 Resor Dray. = 
D Decrease 
ar increase 


f/f. NOWzetien 
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Routing or Preprogram- 
ed Part of Controller 


Array 
wae >, LAV.. 
ee | PBZ 
Complete 
ETS (I)=ETS (I+1) 
No 
/ / 
Call Scan Array [=1,10 


ETC1 (1) 


ETS (10)=ETS TAV (I) 
TEZa) 


Calculate Derivative 









TAV, TB2 
oo 
ee eee 
Teoeneuling 
Logic Tables Part 
\ 
Reactivity Vec oo 
2 
* No 
Change 
es 
2 
No 


Fig. 30. Small deviation controller Dart OL 
digital controller. 
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APPENDIX H 


A. DIGITAL COMPUTER PROGRAM 

The following pages are a copy of the digital computer 
program which was used for the complete Hybrid computer 
Simulation and controller implementation. The program includes 


all parts described in the earlier sections and appendixes. 
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